Electromagnetic structure of A— 2 and 3 nuclei and the nuclear current operator 
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Different modeis for conserved two- and tliree-body eiectromagnetic currents are constructed from 
two- and three-nucleon interactions, using eitlier meson-excliange medianisms or minimal substitu- 
tion in tlie momentum dependence of these interactions. The connection between these two different 
schemes is elucidated. A number of low-energy electronuclear observables, including (i) np radiative 
capture at thermal neutron energies and deuteron photodisintegration at low energies, (ii) nd and 
pd radiative capture reactions, and (iii) isoscalar and isovector magnetic form factors of "^H and ^He, 
are calculated in order to make a comparative study of these models for the current operator. The 

■ realistic Argonne «i8 two-nucleon and Urbana IX or Tucson-Melbourne three-nucleon interactions 
D ' are taken as a case study. For A=3 processes, the bound and continuum wave functions, both below 

pL^ \ and above deuteron breakup threshold, are obtained with the correlated hyperspherical-harmonics 

method. Three-body currents give small but significant contributions to some of the polarization 
observables in the ^H(p, 7)^He process and the2H(n,7)^H cross section at thermal neutron energies. 
It is shown that the use of a current which did not exactly satisfy current conservation with the 
two- and three-nucleon interactions in the Hamiltonian was responsible for some of the discrepancies 
^ ' reported in previous studies between the experimental and theoretical polarization observables in 

QQ ' pd radiative capture. 

o 

CN ■ I. INTRODUCTION 

o ■ 

The present study investigates a number of different models for the nuclear electromagnetic current derived from 
■~»>^ realistic interactions. The emphasis is on constructing two- and three-body currents which satisfy the current conser- 
I vation relation (CCR) with the corresponding two- and three-nucleon interactions. Two different methods are adopted 
' > to achieve this goal: one is based on meson-exchange mechanisms, the other uses minimal substitution in the explicit 
' and implicit-through the isospin exchange operator-momentum dependence of the interactions. A by-product of this 
^ . analysis is, in particular, the elucidation of the sense in which these two different methods are related to each other. 
C ' A variety of electromagnetic observables involving the A=2 and 3 nuclei are taken as case study for these current 
^ , operator models, including the np radiative capture at thermal neutron energy, the deuteron photodisintegration at 
r • low energy, the magnetic form factors of "^He and "^H, and the nd and pd radiative captures. These processes have been 
rN \ extensively studied in the past by several research groups (for a review, see Ref. [1]). Most recently, the authors of 

■ the present paper (and collaborators) have investigated the A—2> radiative capture reactions below deuteron breakup 
" " ' threshold in Refs. [2, 3], and the trinucleon form factors in Ref. [4]. Below, we briefly review those aspects of these 

earlier works that are more pertinent to the present study. 

The ^=3 bound- and scattering-state wave functions were obtained using the pair-correlated hyperspherical har- 
monics (PHH) method [5-7] from a realistic Hamiltonian model consisting of the Argonne uig two-nucleon [8] and 
Urbana IX three-nucleon interaction [9] (AV18/UIX). This technique allows for the inclusion of the Coulomb interac- 
tion in both the bound and scattering states. The nuclear electromagnetic current operator included, in addition to 
one-body convection and spin-magnetization terms, two- and three-body terms. The dominant two-body terms were 
constructed using meson-exchange mechanisms [10] from the momentum-independent part of the AVIS, including the 
long-range pion-exchange component, and coincide with those derived here within the same approach. They satisfy 
the CCR with this part of the interaction. 

The two-body currents originating from the spin-orbit components of the AV18 were constructed using again 
meson-exchange mechanisms [11] [a and lo exchanges for the isospin- independent, and p exchange for the isospin 
dependent terms), while those from the quadratic momentum-dependent components were obtained by gauging only 
the momentum operators [12], but ignoring the implicit momentum dependence which comes through the isospin 
exchange operator. The resulting currents are not strictly conserved. This lack of current conservation was pointed 
out in Ref. [12], but has not been sufficiently emphasized in subsequent papers, mostly because of the short-range 
character of these currents and their generally small associated contributions to photo- and electronuclear processes. 



2 



for example see Rcfs. [2, 11-13]. The overcome of this hmitation is one of the main aims of this work. 

Earher studies, as well as the present one, also take into account the two-body currents, associated with the p-n^ 
and LOTT^ transition mechanisms and with the excitation of intermediate A resonances (for a review, see again Rcf. [1]). 
However, these currents are purely transverse, and therefore unconstrained by the CCR. They are not the focus of 
the present work. 

The effects of A-isobar degrees of freedom in nuclear clectrowcak processes were studied more thoroughly in Refs. [4, 
14], using two different approximations. One was based on first-order perturbation theory (referred to above), while the 
other retained explicit one- and two-A admixtures in the nuclear wave functions via the transition-correlation-operator 
(TCO) method [14]. This latter approach is inherently non-pcrturbative. In particular, it generates three-body 
currents [4], which are strictly not consistent with the three- nucleon interaction in the Hamiltonian. In the present 
work these currents will be derived directly from the three-nucleon interaction, and will satisfy by construction the 
CCR with it. 

The derived new models for the electromagnetic current arc tested in this paper on a variety of A— 2 and 3 processes. 
The predictions for the np radiative capture and deuteron electrodisintegration cross sections at low energies remain 
practically unchanged and are in agreement with the experimental data. For ^4=3 the situation is more interesting 
since the two-body currents play a very important role. For example, in Refs. [2, 3] it was found that two-body 
currents play a crucial role in reproducing the cross section and polarization observables measured in pd and nd 
radiative captures. However, some significant discrepancies remained unresolved. In the nd case, the theoretical 
prediction for the total cross section at thermal energies exceeds the experimental value by 14%. With the present 
model of the current, the overprediction is reduced to 9%. The origin of this overestimate is still puzzling, particularly 
in view of the fact that the astrophysical S*- factor for the pd radiative capture at zero energy is calculated to be within 
1% [15] of that extrapolated from cross section measurements in the range ~ 2 20 keV [16]. 

In the pd case, the calculated tensor observables T20 and T21 at centcr-of-mass (cm.) energy of 2 MeV [3] were 
found to be at variance with data. In that same work, it was also shown that these observables arc sensitive to 
the small (suppressed) contributions arising from electric dipole transitions between the initial pd P-wiwc scattering 
states with spin channel S = 3/2 and the final ^He bound state. When these contributions were calculated in the 
long-wavelength-approximation (LWA) using the Siegert form of the Ei operator, the resulting tensor observables 
were much closer to the experimental values. Since the year 2000, more accurate PHH wave functions have become 
available for the ^=3 nuclei, and the calculations for pd radiative capture could be extended at energies above deuteron 
breakup threshold [17-19]. In preliminary calculations [17], we found that also at 3.33 MeV the theory could not 
reproduce the precise data for the tensor polarization observables T20 and T21 [20]. In the present work, it will be 
shown how the use of a conserved current indeed removes the discrepancy between theory and experiment for these 
observables. Furthermore, the calculation has been extended up to 20 MeV. 

In Rcf. [4] it was shown that the theoretical predictions for the magnetic form factors of "^Hc and "^H were in 
satisfactory agreement with experimental data at low and moderate values of the momentum transfer. The first 
diffraction region, however, was poorly reproduced by the theoretical calculation, especially in the ^He case. The three- 
body current operators, constructed within the TCO approach, gave only very small contributions. This discrepancy 
is not resolved in the present study. 

Alternative descriptions of the A=2 and 3 electromagnetic processes have been also recently reported. A conserved 
current model was developed by Arenohovel and collaborators [21] and applied to A=2 reactions [22, 23]. Several 
groups are studying electromagnetic processes in the three-nucleon system. In Refs. [24, 25], the nucleons are taken 
as interacting via two- and three-nucleon potentials. The electromagnetic currents are then constructed using the 
meson exchange scheme for satisfying the CCR, but only with a part of the interaction. In Ref. [26], the meson 
exchange currents are taken into account using the Siegert's theorem [27]. No three-body currents are considered in 
these works. In Ref. [28], a nuclear model is employed which allows for the excitation of a nucleon to a A isobar, 
and the two-body forces and currents are generated by the exchange of mesons. The A excitation yields also effective 
three-body forces and three-body currents. However, this current model does not satisfy exactly the CCR with the 
adopted Hamiltonian as discussed in Ref. [28]. Very recently, models of the currents derived from chiral Lagrangians 
are starting to appear [29]. In all these calculations, the Coulomb interaction between protons in the scattering 
state is disregarded. Note that, in spite of the differences of the various approaches, the theoretical predictions of 
Refs. [25, 26, 28] and of this paper turn out to be, for most of the observables, quantitatively quite similar. An 
example will be presented for nd radiative capture, where the theoretical results are free from the uncertainty related 
to the omission of the Coulomb interaction. Also other approaches, as the Lorentz integral transform technique [30], 
have been applied to study the electromagnetic response of the trinucleon systems. 

This paper is organized into six sections and three appendices. In Sec. II and III we discuss the model for the 
two- and three-body current operators, respectively. In Sec. IV, wc briefly review the PHH method for the pd and 
nd scattering-state wave function, below and above deuteron breakup threshold. In Sec. V we present results for 
the magnetic structure of the A=3 nuclei, for the deuteron photodisintegration cross section at low energy, for the 
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np radiative capture at thermal neutron energies, and for nd and pd radiative captures at cm. energies up to 20 
MeV. Finally, in Sec. VI, we summarize our conclusions. The connection between the meson-exchange and minimal- 
substitution approaches is elaborated in Appendix A, while a collection of formulas for the two-body current operators 
associated with the quadratic momentum-dependent terms of the two-nucleon interaction, and for the three-body 
current operators in configuration space, are given in the Appendices B and C. 

II. TWO-BODY CURRENT 

The nuclear electromagnetic charge, p{q), and current, j{q), operators can be written as sums of one-, two-, and 
many-body terms that operate on the nucleon degrees of freedom: 



p{q) = E '^^(«) + E P'^^ (9) + • • • , (2-1) 

i i<j 

3{q) = Y^jiiQ) + Y^iiAq) + • • • • (2.2) 

i i<j 

The one-body operators pi (q) and ji (q) are derived from the non-relativistic reduction of the covariant single- nucleon 

current, by expanding in powers of 1/m, m being the nucleon mass. In the notation of Ref. [1], the one-body charge 
operator in configuration space is given by 

Pi{q) = Pi,NR(9) + pt,Rc{q) , (2.3) 
where the leading order term, labeled NR, is 

Pi,Mq) = ^i^'"'"' . (2-4) 

with 

^i = \[G%{ql) + Gl{ql)T,,.] , (2.5) 

while the term labeled RC is proportional to and is explicitly listed in Ref. [1]. In Eq. (2.5) (q^) and G^{q^) 

are the isoscalar and isovector combinations of the nucleon electric Sachs form factors, respectively, evaluated at the 
four- momentum transfer g^. 

The electromagnetic current operator must satisfy the current conservation relation (CCR) 

q-j{q) = [H,p{q)], (2.6) 

where the nuclear Hamiltonian H is taken to consist, quite generally, of two- and three-body interactions, denoted as 
Vij and Vijk respectively, 

^ = E||+E-^.+ E y^^^- (2-7) 

i i<j i<j<k 

Realistic models for these interactions contain isospin- and momentum-dependent terms which do not commute with 
the charge operators. To lowest order in 1/m, Eq. (2.6) separates into 



q-ji{q) = 



|^,P,NR(g) 



(2.8) 



Q-3ij{l) = hj^Ph^niQ) + Pj,NR{q)] , (2.9) 

and similarly for the three-body current jijk{q)- It has been tacitly assumed that two-body terms in p{q) are of order 
Xjrr?. The one-body current is easily shown to satisfy Eq. (2.8). However, it is rather difficult to construct conserved 
two- and three-body currents. 
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It is useful to adopt the classification scheme of Ref. [31], and separate the current jij{q) into model-independent 
(MI) and model-dependent (MD) parts, 



jii(g)=jT(g)+4'^(g)- (2.10) 

The MI two-body current Jij^ {q) has a longitudinal component, constructed so as to satisfy the CCR of Eq. (2.9) (see 
following subsections), while the MD two-body current j^^{q) is purely transverse and therefore is un-constrained 
by the CCR. The latter will not be discussed any further in the present section; it suffices to say that it is taken to 
consist of the isoscalar p7r7 and isovector ujt^^ transition currents, as well as of the isovector current associated with 
excitation of intermediate A resonances [2-4] . 

A method to derive Jif^ {q) was developed by Riska and collaborators [10, 12, 32, 33] and Arenhovel and collab- 
orators [21] (see Ref. [1] for a review). An alternative approach, which we will revisit and generalize in the present 
work, is based on ideas first proposed by Sachs in Ref. [34], and later applied by Nyman in Ref. [35] to derive the 
magnetic-dipole transition operator due to the one-pion-exchange potential. We will refer to these two different ap- 
proaches as the "meson-exchange" (ME) and "minimal-substitution" (MS) scheme, respectively. To appreciate the 
differences and similarities between them, they are discussed in the two following subsections. 

In the rest of the paper, we will use the following notation: a generic nucleon-nucleon interaction will be written as 



^^3=^l^+<f, vlf = <j+<,. (2.11) 

where vf^ and vl^ are the isospin-symmetry conserving {IC) and breaking {IB) parts of the potential, respectively. 
In turn, and are the momentum-independent and momentum-dependent components oivlj'' , respectively. The 
next two subsections deal with the part of the potential. The two-body current associated with its vj^ part will 
be considered in Sec. II C. 

For later reference in Sec. V, a summary of the models for two-body current operators used in the present work is 
given in Sec. II D. 

A. The two-body current operator in the meson-exchange scheme 
First consider the isospin-conserving momentum-independent part of the potential , which can be written as 

""ij = ■^l.ij + ^2,ij Ti ■ Tj , (2.12) 

where Tj and Tj are the isospin Pauli matrices, and vi and V2 are in general functions of the positions and spin operators 
of the two nucleons; V2 includes the long-range one-pion-exchange component. In particular, the isospin-dependent 
terms are given by 

V2.i] n ■ Tj = [Vr{ri;j) + V„r{rij)(T^ ■ (Tj + Vtrinj) Sij] Ti ■ Tj , (2.13) 

where Sij is the standard tensor operator, C7j and aj are the spin Pauli matrices, and the notation of Ref. [36] is used. 
In momentum-space, V2 reads 

V2{k) Ti ■ Tj = [Vr{k) + Var{k) fc^ CTj • aj + Vtr{k)Sij{k)]Ti ■ Tj , (2-14) 

where Vr{k), v„r{k) and Wtr(fc) are related to their configuration-space counter-parts by the relations 

/■OO 

v^{k) = Att I r^dr jo{kr)vr{r) , (2-15) 
Jo 

At: f°° 

''•^^{k) = J r^dr[jo{kr)-l]v„r{r) , (2.16) 
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The factor [jaikr) — 1] in the expression for v^rik) ensures that the volume integral of t'crr(^) vanishes [12], and the 
tensor operator in momentum-space is defined as Sij(k) = k^{cri ■ CTj) — 3(crj • k){crj ■ k). 

If the isospin-dependent interaction V2{k) is assumed to be induced by n- and p-meson exchanges, as for example 
in the Bonn model [37], then 



V2{k) Ti ■ Tj = Vps{k) + [2vp{k) + U7r(fc)] k'^ CTi ■ a-j + [vp{k) - VTr{k)]Sij{k) 



Ti ■ Tj , 



with 



Vnik) 

Vp{k) 

VpS 



J-kNN 



Sm'i k"^ + mi 



91nn{'^ + i^pNN? fp{k) 



I2rn? 



k'^ + m?p ' 



(2.18) 



(2.19) 
(2.20) 
(2.21) 



where fjrNN, and rjpNN and K.pNN are the coupling constants of the tt- and /?- mesons, f-,t{k) and ,fp{k) are the associated 
form factors (usually, of monopole type), and nip are their masses, and finally m is the nucleon mass. 

More generally, if one assumes that the interaction V2{k) is due to the exchange of a number of "7r-like" pseudoscalar 
{PS) and "/9-like" vector (V) mesons, then one finds 



V2{k) Ti ■ Tj = vvs{k) + [2vv{k) + vps{k)\ fc^ cTj • (Tj + [vv{k) - vps{k)]Sij{k) 
where the functions vps{k), vv{k), and vvs{k) are given by 



n ■ Ti 



(2.22) 



N 



N 



vv{k) = 
vvs{k) = £/ys,ap^ 



+ 



a=l 

JV 



a=l 



(2.23) 
(2.24) 
(2.25) 



In the expressions above, rua is the mass and / 



PS,a 



-flNNj^rnl, = -5^ArA,,a(l + ^pN N .af /I2m? and 



fvs,a = 9pNN,a Coupling constants of the exchanged a-meson. These parameters are fixed so that 



vps{k) = [var{k)-2vtr{k)]/i 

Vv{k) = [v„r{k) + Vtr{k)]/3 , 
vvs{k) = Vr{k) , 



(2.26) 
(2.27) 
(2.28) 



where v^r{k), Vtrik), and Vr{k) are given in Eqs. (2.15)-(2.17). 

The two-body currents due to these PS- and T^-mcson exchanges are then derived by minimal substitution in 
the effective PS-NN and V-NN coupling Lagrangians. The non-relativistic reduction of the associated Feynman 
amplitudes in momentum space leads to: 



jij {ki , kj ; PS) = 3i Gl {ql ) (Ti XTj), 



vps{kj)(Ti{aj ■ kj) - vps{ki)crj{ai ■ fcj) 



+ 



kf - k] 



[vpsih) - vps{kj)] {(Ti ■ ki){(Tj ■ kj) 



(2.29) 
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jij {ki,kj;V) = - 3i Gl {ql ) {n x t,- 



vv{kj)(Ti X ((Tj X kj) — vv{ki)<Tj X (cr, x fe,) 



_ f^2 [(^» k,j){ai X fej) • (<Tj X fej) 



+(crj X fei) cTj • [ki X fej) + (cTj X kj) cTi ■ [ki x fcj)] 



(2.30) 

ii,-(fei, fe,; = -iGl{ql){n X T,). §^ [vvs{ki) - vvs{kj)] , (2.31) 



where ki and fej are the fractional momenta dehvered to nucleon i and j, with ki + kj = q, and G^{q'^) is the isovector 
combination of the nucleon electric Sachs form factors [1]. The current 

3ij^ {ki,kj\xP) = jij{ki,kj; PS) + jij[ki,kj;V) -\- jij{ki,kj;VS) , (2.32) 

satisfies exactly the CCR with the potential given in Eq. (2.12). 
Configuration-space expressions are obtained from 

Mq;B) = J dxe^'^-J ^^e^'=.-(^.-)ei'=H-.-)jV fe,; B) , (2.33) 

where B=PS, V or VS, and can be found in the Appendix of Ref. [12]. For later reference, we report below the 
configuration-space expression for the current associated with the isospin-dependent central potential jij{q; VS): 

jij{q;VS) = -Gl{ql){nxrj)^e^'^-^f2^-^ /v/^^ + ^) '"''"'^ '"'^"^'^ ' 

L^ix) = ^ml + ^(l-4x^) , (2.34) 

where r = — rj and R= {ri + rj)/2. 

The constriiction of the two-body currents associated with the isospin-conscrving momentum-dependent part of 
the interaction v'^j is less straightforward. A procedure similar to the one just reviewed has been applied to the 
case of the currents from the spin-orbit components of the interaction [11]. It consists, in essence, of attributing 
these to exchanges of "cr-like" and "w-like" mesons for the isospin-independent terms, and to "/o-like" mesons for the 
isospin-dependent ones. Explicit expressions for the resulting currents can be found in Ref. [11]. 

The two-body currents from the quadratic momentum-dependent terms of the interaction are listed in Ref. [12] 
and were obtained by minimal substitution, that is Pi ^ Pi — ei A{ri). While minimal substitution ensures current 
conservation for the isospin-independent (quadratic-momentum-dependent) components of the interaction, this pre- 
scription does not lead to a conserved current for the isospin-dependent ones. Indeed, the commutator in Eq. (2.9) 
gives rise to an isovector term proportional to (xj x Tj)z, which cannot be generated by minimal substitution (for a 
discussion of this point, see Ref. [38]). These isovector currents were ignored in all previous work since, in view of 
their short range, they were expected to give negligible contributions. 

The currents associated with vl^ , which have never been considered up until now, will be discussed in Sec. II C. 

B. The two-body current operator in the minimal-substitution scheme 

Consider again the isospin-conserving momentum-independent part of the potential given in Eqs. (2.12) and (2.13). 
The isospin operator Tj • Tj is formally equivalent to an implicit momentum dependence [34] , since it can be expressed 
in terms of the space-exchange operator, P^- , using the formula 

Ti-Tj = -l-{l + <Ji- <jj)Pij , (2.35) 

valid when operating on antisymmetric wave functions, as in the case of a fermionic system. The space-exchange 
operator is defined as 
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P,,- fin, rj) = e^oi-^^+nr^^ f{rurj) = /{rj, n) , (2.36) 

where the V-operators act only on the generic function f{ri,rj) and not on the vectors = — rj = —rji in 
the exponential. In the presence of an electromagnetic field, after performing minimal substitution, the operator 

becomes [34] 



P 



= q'^H ■ ['Vi-i€iA{ri )]+rij ■ [Vj-ie-iA{rj )] 



(2.37) 



where A{r) is the vector potential, and the functions gi{ri) and gjivj) have been defined as gi{ri) = —icirji ■ A{ri) 
and gj{rj) = —iej rij ■ A{rj). The operator P-j is then the product of two operators, each having the general form 



P(r) 



„a-V+g(T-) 



where a is a vector independent of r. It has been shown in Ref. [34] that P(r) can expressed as 



P(r) 



go Jr 



(2.38) 



(2.39) 



where ds is an infinitesimal element of a straight line parallel to a, which goes from position r to position r + a. Using 
this general result in Eq. (2.37), we obtain 



■^ij — e 3 Jr'ij , 



(2.40) 



with ds = ds Vji and ds' = ds Vij . The line integrals arc performed on straight lines parallel to rji and f ^ . 

The procedure of Ref. [34] leading to Eq. (2.40) can be generalized and the two integrals can be performed on two 
generic paths 7ij and 7^^, that go from position i to position j and vice versa, as shown in Fig. 1. Indeed, for a gauge 
transformation 



A{r) -> A{r) - VG(r) , 



(2.41) 
(2.42) 



where G(r) is a generic function, it can be shown that the state 



l_(l + ^..^.)p. 



A 

3)^ ij 



where P^ is given in Eq. (2.40) with the generic integration paths of Fig. 1, transforms as 



(2.43) 



l-(l + o-,-(T,)i^f 



l-(l+0-,-<T,)P, 



A-VG 

J)-" ij 



l-{l + cri-aj)P,^ 



(2.44) 



in conformity with the requirements of gauge invariance of the theory. A conserved current can then be derived by 
considering an infinitesimal gauge transformation, as discussed in Ref. [34]. 

Rather than following the general procedure of Ref. [34] , we obtain the current in the limit of weak electromagnetic 
fields, since calculations of photo- and electronuclear observables are typically carried out in first-order perturbation 
theory in these fields. Then, by retaining only linear terms in the vector potential, we find 
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— ici / ds ■ A{s) — icj / ds' ■ A{s') 
= v% - j j,j{x) ■ A{x)dx , 



^ Vij+V2,ij 



(2.45) 



where the paths •jij and are those illustrated in Fig. 1, and in the second and third lines of the equation above v^j 
is as defined in Eq. (2.12). The current density operator is then given by 



jij{x) = iV2,ij 

and its Fourier transform reads 



Cj / dsS{x — s) + Cj / ds' S{x — s') 



{1 + Ti-Tj) , 



(2.46) 



jij (q) = i V2,ij ( / ds e"''^ + ej / ds' ( 



(1 + Ti-T,). (2.47) 
A number of comments are now in order. Firstly, the current in Eq. (2.47) by construction satisfies the CCR 



q-jiAo) = 



V-',-,Pi,NR(q) +Pj,NR(g) 



(2.48) 



This is easily verified by observing that along any path from to Vj 



and that 



ig- / ' ds e"i-^ = e"^ ''' - e"i ''' 



i {ei - ej){l + Ti ■ Tj) = GM) in X Tj) 



(2.49) 



(2.50) 



Secondly, because of the arbitrariness of the two integration paths, the prescription just outlined does not lead to 
a unique two-body current. Moreover, if V2,ij consists of a sum of different terms, then, for each of these, different 
paths Jij and 7^-^ can be selected. Hence, Eq. (2.47) can be interpreted as a parameterization of all possible two-body 
currents which satisfy the CCR with the two-body potential given in Eq. (2.12). In particular, it is interesting to 
note that the longitudinal part of the two-body currents of Eqs. (2.29)-(2.34) obtained in the ME approach can also 
be derived in the MS scheme. This connection is shown in Appendix A. 

Thirdly, the choice of a linear path {LP) for -y^j and 7^^ in Eq. (2.47) (as in Ref. [34]), namely 



s = —s 



xr. 



I] 1 



< a; < 1 , 



gives 



with 



Ji^il) = '^V2,ij Gl{ql) in X Tj)z Tij fij{q) , 



(2.51) 



(2.52) 



fijiq) = 



q ■ Tij 



(2.53) 
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Note that /„ (q = 0) = i. 

Lastly, it is important to observe that, in the hmit q ^ 0, the current operator jij{q) becomes path-independent, 
i.e. it is unique, and is given by 



jij{q = 0) = -V2,ijGE{qi,) {n x Tj), nj . 
This result can also be derived in a more direct way by considering the following identities 



(2.54) 



jij{q = 0) = j dxjij{x) = - J dxxV ■ jij{x) 



(2.55) 



where in the first line the volume integral of jij{x) has been re-expressed in terms of the divergence of the current, 
ignoring vanishing surface contributions, and in the second line use has been made of the CCR. Evaluation of the 
commutator leads to Eq. (2.54) above. Note that, mutatis mutandis, namely jij{x) j{x) and V2,ij H, etc., 
the second line of Eq. (2.55) is, in essence, the Siegert theorem for the electric dipole operator, to which j{q = 0) is 
proportional. 

We now derive the current operators associated with the momentum-dependent operators of the two-nucleon inter- 
action, within the present scheme. In the case of the spin-orbit interactions, and V2 of Eq. (2.12) are 



Vb{r) L - S , 
VbT{r)L- S , 



(2.56) 



where the notation of Rcfs. [8, 36] is used, and L = Tij x (pi — Pj)/2, pi and pj being the particles' momentum 
operators. Performing minimal substitution in vi, we obtain 



jij{q; h) = -Vb{r) ( ei e^"'^' - ej e^" "^^ ] S x r,, . 



(2.57) 



For the isospin-dependent term uf, we first symmetrize as 



""tij '^i ■ "^J = 9^'6r(r-) {L- STi- Tj + Ti ■ Tj L ■ S) 



(2.58) 



and then perform minimal substitution in both L ■ S and r, • Tj . When only terms linear in the vector potential A 
are kept and the path 7^^ = —-jij is taken, the associated current is found to be 



3ij{q;hT) = -Vbr{r)S x Tij ( rjje''^''' - rjie"^''' 



+ :;Vbr{r)Gl{q^^){Ti X Tj)4 L S / dse"^-^+ / dse'"'^ L-S ] , 



/ 

with r]i = G§(9^) Ti ■ Tj + G\{(^^Ti^z- In particular, the linear path of Eq. (2.51) leads to 



(2.59) 



j5-^(9; br) = \vbAT)S X Vij ( r^.-e'"-* - ^e'"- 



+ '7)'VbT{r)G\{ql)(Ti X Tj), 



L ■ S Tij fij (q) + Tij fij {q)L- S 



(2.60) 



with fij{q) defined in Eq. (2.53). 

The current operators associated with the quadratic momentum-dependent terms of the interaction can be derived 
in a similar fashion. Their explicit expressions are listed in Appendix B. 
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C. Two-body current associated with the isospin-symmetry breaking interactions 

The current operators constructed so far in the ME and MS schemes satisfy the CCR with the isospin symmetric 

component of the two-nuclcon potential. However, the latest realistic models of the nucleon-nucleon interaction 
contain isospin-symmetry breaking terms. In the notation of Ref. [8], this part is written as 

< = E ^pi^i^n, , (2.61) 

p=15 

and the four isospin-symmetry breaking operators have the form: 

QP-15,...,18 ^ j,^^ ^ ^, _ ^^y,^ ^ rp^^ ^ _l_ ^ (2.62) 

where Sij is the standard tensor operator and the isotensor operator is defined as = 3Ti.2 Tj.^ — • Tj. The 
dependence on Xj • Tj generates two-body currents which can be taken into account by modifying the isospin-dependent 
central, spin-spin and tensor terms of the potential as 



Vr{r) = [vr{r) - vx5{r)] , 

Varir) = [Var{r) - ViQ{r)] , 

vtr{r) = [vtrir) - virir)] , (2.63) 

and by using Vrir), v^rir), and Vtrif), instead of Vr{r), v„r{f), and Vtrir) in V2,ij of Eq. (2.13). However, the con- 
tributions associated with the currents from these isospin-symmetry breaking terms have been found to be negligibly 
small for all the observables of interest here. 



D. Summctry of two-body current models 



For the sake of clarity and for later reference, we summarize here the salient features of the three different models 
for the model-independent current jf^^ {q) corresponding to the AV18 interaction [8], considered in the present paper. 

• Old-ME model. This model is that introduced in Refs. [2-4] and discussed in Sec. II A. It is given by 



= v') + j^^iq; vn . (2.64) 

We re-emphasize that, while jij^{q;v^) satisfies the CCR with (which includes the long-range one-pion- 
exchange term), 3if^{q;v'P) is not strictly conserved, as discussed in Sec. II A. 

• New-ME model. This model retains 3if^{q; u°) for the momentum-independent interaction, as in the "old-ME" 
model. For the momentum-dependent interaction, it uses instead the two-body current obtained in the MS 
scheme with a linear path, explicitly 



^.M/,new(^) ^^.M^^^^.^O) +_,;^P(q.„P) , (2.65) 

where 



3ii{q; = Jio (g; b) + j;§^(g; br) + jij {q; LL) + j;§^(q; LLr) 

+ jij{q\bh)+3^^{q;bbT) , (2.66) 

and jij{q;b), jtj^{q;bT), j,j{q:LL), j[/{q]LLT), jtj{q;bb) and j^^{q;bbT) are listed respectively in 
Eqs. (2.57), (2.60), (B5), (B6), (BIO), (Bll). In addition, the isospin-symmetry-breaking contributions are 
included via Eq. (2.63). The two-body current operator given above satisfies exactly the CCR with the AV18 
potential. 

It is important to stress that the longitudinal component of j^^{q; v^) can also be obtained in the MS scheme, 
as discussed in the previous section and in Appendix A. 
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• Linear path MS (LP-MS) model. This model uses a two-body current obtained in the MS scheme using a linear 
path, exphcitly 

= j-r (g; v') + j-r (9; vn , (2.67) 

where jij^{q;v°) is the current given in Eq. (2.52) and jij^{q;vP) is the same as in Eq. (2.66). 

III. THREE-BODY CURRENT 

Thrcc-body currents involving the excitation of an intermediate A resonance were derived recently in the context 
of a study of explicit A components in the trinucleon wave functions [4, 28]. In addition, the three-body current 
associated with 5- wave wN scattering on an intermediate nucleon was also included in Ref. [4]. The conclusion of 
that work was that these tlirec-bocly mechanisms give a small contribution to the magnetic form factors of "'ll and ''He 
over a wide range of momentum transfer. However, the three-body currents considered in Ref. [4] were not strictly 
consistent with the three-nucleon interaction (TNI) included in the Hamiltonian. 

In this section we generalize the meson-exchange (ME) and minimal substitution (MS) approaches to the case of 
the three-body current induced by a TNI Vijk- The resulting current satisfies, by construction, the CCR with Vijk- 
To be specific, we consider below the Urbana-IX model [9], but the methods that are developed are applicable to 
other phenomenological models of TNI's, such as the Tucson-Melbourne [39] and Brazil [40] ones. 

A. The three-body current in the meson-exchange scheme 

The Urbana-type TNI is written as sum of a short-range spin- and isospin-independent term and a term involving 
the excitation of an intermediate A resonance. The central term is irrelevant to the following discussion and is 
therefore ignored, while the A-excitation term is given by [9] 

^ijk = ^ Vj;ki , (3.1) 
cyclic ijfe 

Vj;ki = ^27r^{^ij i Xjf.} {Xj ■ Tj , Tj ■ Tfe} -|- - [Xij , Xjk] [Xj ■ Tj , Tj ■ Tfc] ^ , (3.2) 

where {. . .} ([. . .]) denotes the anticommutator (commutator), 



Xij = vii{r)ai ■ a, + <^(r)5,, , (3.3) 

and v^^^{r) and «//(?') are the standard spin-isospin and tensor-isospin functions occurring in the one-pion-exchange 
interaction, modified by a short-range gaussian cutoff. The parameter ^2,^ as well as the strength of the central 
term are determined by reproducing the triton binding energy in a Green's function Monte Carlo calculation, and the 
nuclear matter equilibrium density in an approximate hypernetted-chain variational calculation [9]. 
In momentum space, Vj-ki can be conveniently expressed as 



6{ki + kj + kk) , (3.4) 
where the A^A-transition interaction is defined as 



Vj-ki{kj, ki, kk) — — 



u]fc(fefc; 



AN NN) Vij {ki; NN NA) + h. c. 



Uij{k;NN ^ NA) 



.II 



{k)S^I{k) 



(3.5) 



Here S. 



and Tj are the spin- and isospin-transition operators that convert nucleon j into a A-isobar, and S-j{k) is 



the momentum-space tensor operator in which the Pauli spin operator of particle j is replaced by S, . The functions 



v^i^{k) and vl^{k) are related to their configuration-space counter-parts by similar relations to those in Eqs. (2.16) 



12 



and (2.17). The momentum transfers to nucleons i, j, fc, respectively ki, kj, and kk, sum up to zero. Manipulation 
of products of transition spin and/or isospin operators is facilitated by making use of the following identity: 

S/3 = -Sa0 - gCa/S^cr^ . (3.6) 

The A'^A-transition interaction is assumed to originate from exchanges of "7r-like" (PS) and "p-likc" (V) mesons, 
with the associated components Vpg{k) and Vy (k) related to v^^^{k) and vj.^{k) by relations identical to those in 
Eqs. (2.26) and (2.27). Thus, the PS- and F-exchange three-body currents, illustrated in Fig. 2, in momentum space 
read 



if,f(q) = I A 



2 



2tt 



v^ikk-, AN ^ NN) \ji;{ku kj; PS) + jlj (k^ kj; V)] 



+ [jH {kj , fcfc ; PS) + j]l (fc, , kk ; V)] \ij {ki ; NN ^ iVA) + h. c. 



(3.7) 



where the PS and V currents jl^ {ki, kj: PS) and jlj {ki. kj; V) involving A excitation are obtained from Eqs. (2.29) 
and (2.30) with the replacements aj {xj) — + Sj {Tj). Configuration-space expressions are listed in Appendix C. 
Finally, the current above satisfies the CCR with the TNI of Eq. (3.4). 

B. The three-body current in the minimal-substitution scheme 

Consider the isospin dependence of the TNI. The anticommutator term is first expressed as 

{t, ■ Tj , Tj • Tfe} = 2 Tj • Tfe , (3.8) 

and the associated current j^^iil) ii^Gu derived with the same methods discussed in Sec. II B and is given by (see 
Eq. (2.47)) 



j^^M = 2i^2. {Xij , X,k} U [ dse"^-' + Ck [ ds' e"^ A {1 + n- tu 



(3.9) 



where jik and are generic paths from to rfe and rk to (see Fig. 3). 
In the case of the commutator term, we first note that 

\ [Ti ■ Tj , Tj ■ Tfe] = PljPJk - PJkPlj , (3.10) 

where P-^ = (l + ■ Tj)/2 is the isospin-exchange operator. The product PljPJ^, when acting on antisymmetric wave 
functions, is equivalent to 

PijPjk = PjkPijPjkPij 1 (3-11) 

where P^j is the spin-exchange operator, defined similarly as P^^-, and is space-exchange operator introduced in 
Eq. (2.36). Note the ordering of the operators on the r.h.s. of the equation above. Obviously, the product PJkPij 
is given by a relation similar to Eq. (3.11) in which the order of the jk and ij pairs is inverted. The products of 
space exchange operators PjkPij and PijPjk are equivalent, respectively, to the exchanges {ri,rj,rk) {rk,ri,rj) 
and (rj, Vj, rk) ^ {fj,rk, Vi) (see Fig. 4), and can formally be expressed by the operators 

PjkPij = e^*''-'^*+''«-'^i+^i'=-'^'= , (3.12) 
PijPjk = e^i*-'^'+^*'^-'^^ , (3.13) 
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where, as before, the gradients do not act on the position coordinates in the exponential. The methods of Sec. II B 
can now be apphed to the present case. Gauging the gradient operators and retaining only linear terms in the vector 
potential (valid for weak electromagnetic fields) lead to the following current jf-hiiq) from the commutator term of 
the TNI: 



jf-kiio) = -^MA^ij , Xjk] 



€i I dse'i-^ + Cj [ dse'i-^ + Ckf ds c'" '' ){1 + n ■ Tj){l + tj ■ Tk) 



( 

- (a I dse'^-' + Cj I dse'i-' + ekf ds e'^'^^il + Tj ■ Tk){l + n ■ tj) 

^ J0U J^'ik ' 

where /3jfc (/3fei) is a generic path starting at (r^) and ending at (r^), and so on (see Fig. 3). 
The expressions for j^j^. ■ (q) and j^j.^ (q) may be simplified by selecting the following paths: 

7fei = -lik , P'ij = -fiji , (i'jk = -Pkj , P'ki = -f^ik ) Pik = -Pji - Pkj , 



(3.14) 



(3.15) 



namely the path from rk to r, is taken to be the same as that from n to rk but in the opposite direction, and so 

on. The last relation means that the path f3ik from to is chosen to go through the position rj exactly along the 
paths —Pji and —Pkj (the two latter paths are still arbitrary). We then obtain: 



= 2 Gl{ql) {Xi, , X,k} in X n). / 



dse" 



{Ti,z Tj ■ Tk - Tj^z Ti ■ Tk) / ds e"^'^ 
+ (Tk,z Ti ■ Tj - Tj^z n ■ Tk) / ds e"^'^ 



L 



(3.16) 



This current is easily shown to satisfy the CCR with the TNI. As in the case of two-body currents, the limit q=0 
is path-independent. 



3j;kiil = 0) = i 



j;ki 



J dxx [pi,NR{x) + pj,NR{x) + pk,NR{x)] 



(3.17) 



Furthermore, when the paths ^ik, Pjk and /3^ - are taken as straight lines as in Eq.(2.51), then Eq. (3.16) becomes 



{Ti,z Tj ■ Tk -Tj^z Ti ■ Tk)rjk fjk{q) 
+ {Tk,z Ti ■ Tj - Tj^z Ti ■ TkjVij fij{q) 



(3.18) 



where the functions fij{q) are defined in Eq. (2.53). 

Finally, note that the present approach can also be used to derive the currents associated to the the Tucson- 
Melbourne (TM) [39] and Brazil [40] TNI interaction models, since these can be cast in the form [41] 



Vj-ki = Fs{j; ki) {Ti ■ Tj , Tj ■ Tk} + FaU; ki) [n ■ Tj , Tj ■ Tk] 
For example, the TM model has 



(3.19) 
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Fs{j; ki) = cs {X,j , Xjk} + B{rij,rjk) {Sij + <Ji ■ <jj , Sjk +CTj-<Jk) , (3.20) 

FA{j;ki) = ca[X,j,Xj,] , (3.21) 

where the parameters cs and ca have the values cs — —0.063 and ca — —0.018 [39], and the function B{rij,rjk) 
depends on a cutoff A, fitted to reproduce the triton binding energy. 



C. Summary of three-body current models 

Wc summarize in the present subsection tlie different models for the three-body current used in the present study. 

• Old-TCO model. The model is that introduced in Ref. [14] and subsequently refined in Ref. [4]. As already 
mentioned, it does not satisfy the CCR with the Urbana or Tucson-Melbourne TNIs. 

• ME-model. In the case of the Urbana- type TNI, the three-body current jj^f (q) satisfying the CCR is given by 

the configuration-space expression of Eq. (3.7), which can be derived from Eqs. (CI) and (C2). For the TM-type 
TNI, some difficulties arise, since the second term of the operator Fs{j; ki), proportional to B{rij,rjk), cannot 
be simply related to the exchange of a single 7r-like or p-like meson. Therefore, in this case a hybrid approach 
is used, where the current associated with this last term is treated within the linear-path MS scheme, while the 

rest is obtained within the ME scheme. 

• Linear path MS {LP-MS) model. Within the MS scheme, we select again the linear path of Eq. (2.51) to 
construct the three-body current, as given in Eq. (3.18). 

Note that the current corresponding to the TNI defined in Eq. (3.1) involves a cyclic sum over ijk, i.e. 

jijk{q) = J2 jrMi) ■ (3-22) 

cyclic ijk 

Lastly, it is worth remarking here that, at low values of the momentum transfer, the contributions associated with 
the operators jffjfiq) and j^jkiQ) (as well as j^^''^'^^{q) and jif{q)) are calculated to be essentially the same for 
the observables of interest in the present study. 



IV. WAVE FUNCTIONS 



The trinucleon bound-state and N — d scattcring-statc wave functions are obtained variationally with the pair- 
correlated-hyperspherical-harmonics (PHH) method [6]. Recently, in a series of papers [7, 42, 43], the method has 
been generalized to solve the N — d elastic scattering problem above the deuteron breakup threshold (DBT), thus 
allowing for the study of electromagnetic processes at higher energies than previously treated [2, 3]. For completeness, 
the method will be reviewed briefly and a summary of relevant results obtained for N ^ d scattering observables at 
energies above the DBT will be presented. 

The wave function "^1^2'^' for a TV — d elastic scattering state with an incoming relative orbital angular momentum 
L, channel spin S {S = 1/2,3/2) and total angular momentum JJ^, is written as 



where describes the system in two regions: (i) the "core" region where the three particles are close to each 

other and their mutual interactions are large and (ii) the "breakup" region where the three particles are far from each 
other. The other term \1>^'^'^'^^ describes the system in the 1-1-2 "clusterization" asymptotic region, where intercluster 
nuclear interactions are negligible. The function ^^^-^-^^ (for p — d, as an example) is given by 



LSJJ^ 



JJz 



cyclic ijk L'S' 

SLL'Sss'H]r,{rj,pri) - Sis,L'S'{F)H+{rj,pri) 



(4.2) 
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where (j)d is the deuteron wave function, Xi the spin state of nucleon i, Xi and yi the Jacobi vectors defined, respectively, 
and Ui = (2rj — rj — rk)/V^ = \/4/3ri, and p is the magnitude of the relative momentum between 
deuteron and proton. The functions are defined as 

where and are the regular and irregular Coulomb functions, respectively, and rj is the Sommerfeld parameter. 
Note that for n — scattering 77 = 0, and Fl{0,x)/x and Gl{0,x)/x reduce to the regular and irregular spherical 
Bessel functions. The factor (1 — c-i^r^2L+i been introduced to regularize the function G at the origin, and k is 
taken as a variational parameter. The complex parameters S^g j^,g,{E) are the .S-matrix elements which determine 
phase-shifts and (for coupled channels) mixing angles at the cm. energy E = Tc.m. — B2 where B2 = 2.225 MeV is 
the deuteron binding energy and 



Tc.„^.=pV(2M) , (4.4) 

is the N — d cm. kinetic energy, /i being the N — d reduced mass. The sum over L'S' in Eq. (4.2) is over all values 
compatible with a given J and parity. 

The second term ^q"^" of the trial wave function must describe those configurations of the system where the particles 
are close to each other. For large interparticle separations and energies below the DBT, VE'^'^'' goes to zero, whereas 
for higher energies it must reproduce an outgoing three particle state. In terms of the PHH basis, ^'^'^^ is expanded 
as [6], 

AT,, NK{a) 

^^c' = P-"'^ E E ^<^MZ^K , (4.5) 

a=l K=l 

where p = ^/x^ + yf is the hyperradius. The functions ZaK are antisymmetric under the exchange of any two pairs of 
particles and account for the angle-spin-isospin and hyperangle dependence of channel a, K. The hyperangle is defined 
as cos0j = Xi/p. The index a denotes collectively the spectator i and pair jk orbital and spin angular momenta 
and isospins coupled to produce a state with total angular momentum and parity J^, while the index K specifies 
the order of the Jacobi polynomial in the hyperangle. The values of Nc and Nxict) are increased until the desired 
degree of convergence in the quantity of interest is obtained (see the discussion in Sec. V A). In the PHH approach, a 
correlation factor is included in ZaK in order to better take into account those correlations induced by the repulsion 
of the potential at short distances. In this way, the rate of convergence in the Nk{q) expansion is improved very 
significantly (A^/f (a) < 10 in all cases). 

The functions Ua,K{p) are the hyperradial functions to be determined by the variational procedure, once the 
boundary conditions arc specified. In practice, the functions Ua,K{p) are chosen to be regular at the origin (ua,if(0) = 
0) and to have the following behavior as p ^ 00 

Ua,K{p) ^ , E <0 , 

-»Mp) - EE (^"'^""' ^^'"'^'') . ^ka'K'iE) e^<^^ , E 

a' K' 

where Sj^g^^,j^, are the ^-matrix elements for the process 1 + 2— >1 + 1 + 1. The matrices n and c are defined as 



0^ 

— >0 , 
m 



(4.6) 



(4.7) 



naK,a'K' = Hni / dCl ZIj^ Za'K' , CaK,a'K' = Hni / dQ ZIj^ pVc Za',K' , (4.8) 
p^co J /9->oo J 

where Vc is the Coulomb potential energy and dU = (cos 0i)^(sin (f)i)'^d(f)i dxi dyi. Once the above boundary conditions 
are applied, it has been shown that the Kohn variational principle for scattering states is valid also above the DBT 
(for more details, see Ref. [43]). This principle can therefore be used to compute the matrix elements ^ls L'S'i^) 
and Slg ^,j^,{E) and the functions UaK{p) occurring in the expansion of "^c- This is achieved in practice by making 
the functional 
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[■Sis,L'S'm = Sis,L'AE) - V2,imp{^[li'^\H - E\^[li'') , (4.9) 

stationary with respect to variations in the Slg and UaK- 

Phase-shifts and mixing angles for n — d scattering have been obtained from a realistic Hamiltonian model, and 
have been shown to be in excellent agreement with corresponding Faddeev results [44, 45], thus establishing the high 
accuracy of the PHH expansion for this scattering problem. It is important to emphasize that the PHH scheme 
permits the straightforward inclusion of Coulomb distortion effects in the p — d channel. The PHH results for p — d 
elastic scattering are as accurate as those for n — d scattering. 

For example, various p — d observables at Tc.m. = 6.66 MeV predicted by the AV18/UIX model are shown in Fig. 5 
and are found to be in good agreement with the available experimental data [46, 47]. The large discrepancy observed 
for the Ay and iTn observables (the "Aj^ -puzzle") is connected to a not well understood deficiency of the nuclear 
interaction, most likely of present TNI models. Resolving this "Aj,-puzzle" is a current and important area of research. 

The bound-state wave function "^'^'^^ (,/ = 1/2) is just given by the term '^'^f', which is expanded as in Eq. (4.5). 
In this case, the functions Ua^xip) arc determined by the Rayleigh-Ritz variational principle, applying the boundary 
conditions Ua,K{p oo) 0. In practice, they are expanded in terms of Laguerre polynomials multiplied for an 
exponential factor [5, 6]. The number of channels included in such an expansion will be denoted by iVc(b.s.) in the 
following. The PHH expansion is very accurate also for boiuid states, as shown, for example, in Ref. [48], where a 
very detailed comparison with the results of the Faddeev calculations of the Bochum group has been performed. 

In the following, it is convenient to use the wave function "^'^a2cr, where C72 and a are the spin projections of the d 
and N clusters, and p their relative momentum in the incident channel. It is given by 



1 plCTi 

*il<.=4^E<l'^2,2^l^^.) E ■^'^ {SS.,LM\JJ,)Y£^{p)—^ffi'' , (4.10) 

SS^ LMJJ^ 

where is the Coulomb phase shift. For a nd state the factor e'"^^ is omitted. The wave function ^^l^cr satisfies 
outgoing wave boundary conditions, and is normalized to unit flux, while the two- and three-nucleon bound-state 
wave functions are normalized to one. 

In earlier papers [2, 3], the sum over J in Eq. (4.10) was truncated to a given value Jmax = 7/2, since the analysis 
was limited to study low-energy radiative capture {Tc.m. < 2 MeV). In the present work, we extend the calculations 
to higher energy. In this case, it is necessary to take into account also the contribution of higher partial waves. For 
large values of J, and correspondingly large values of L, the centrifugal barrier between the deuteron and the third 
nucleon prevents the two clusters to approach each other. The corresponding ^'f^/"^" of either the nd or pd state can 
therefore be approximated to describe the free or Coulomb-distorted motion. For pd as an example. 



*f^/"^-*f^/,^;e=4 E \[UXi)®Xi\s'®Yi^'{yi)\^^ ^"^^y^^ J>Jma.. (4.11) 

y o ,. . ^ JJz P'i 

cyclic ijh 

In the calculation of transition matrix elements, we found it convenient to divide the sum over J in Eq. (4.10) as 
J2j J2j<j +X]j>j ■ the first sum, the wave functions Vl/f^/"^' are calculated by taking into account the 
full PHH expansion (the effect of increasing the number A^c(s.s.) of scattering state channels is studied for a few 
selected cases in Sec. VA). In the second sum, the wave function is approximated as in Eq. (4.11). In this 

case, the sum over L can be evaluated analytically to reconstruct the Coulomb distorted "plane wave" describing p — d 
motion. In summary. 



*ila = 4^ E E(^^'1^2|55,)^i^(55,,XM|JJ,)y;M(p)^*[f2''^+A*(-|i,. , (4.12) 

J<J^^^,J^ SS^ LM 



where 



cyclic ij/c J<Jina.x,Jz SSz 

^i^(55„iM|JJ,)y^V(p)^*i+/,L ' (4-13) 



X 

LM 
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and tpc (p, r) is the solution of the three-dimensional Schrodinger equation with the pure Coulomb potential behaving 
asymptotically as a plane plus a scattered wave, i.e. 

ij4+\p,r)=e'P-''e-''''/'^T{l + ir])iFi{-ir],ipr-ip-r) . (4.14) 

Here, iFi is the confluent hypergeometric function. The function given in Eq. (4.14) reduces simply to the plane wave 
for 77 = {nd case). 

V. RESULTS 

In the present section we report results for the np radiative capture at thermal neutron energies and deuteron 
photodisintegration cross section at low energy, the nd and pd radiative capture reactions at cm. energy Tc.m. = 2-20 
MeV, and the isoscalar and isovector magnetic form factors of and ^He. In the next two subsections, we report the 
results for the pd radiative capture at Tc.m. = 2.0 and 3.33 MeV, for which there are very accurate cross-section and 
polarization data [20, 49]. One reason for doing so is to test the quality of the boimd and scattering wave functions, 
in particular by studying the rate of convergence of calculated reduced-matrix elements (RMEs) with respect to the 
number of channels included in the PHH expansions of these wave functions. 

The second reason is to make a comparative study of the different current operator models introduced in the present 
work. Some of these models satisfy the CCR exactly, others do so only approximately. The question is how critical is 
this lack of current conservation and how large are the contributions of three-body currents induced by the trinucleon 
interaction. 

In Ref. [3] significant deviations were obtained between the measured and calculated tensor observables T20 and 
T21 in pd radiative capture. That earlier study was carried out with a current operator including, in addition to 
one-body, two- and three-body terms, denoted as old-ME and old-TCO in the present work. As shown below, most 
of the observed discrepancy between theory and experiment can be traced back to the fact that the current of Ref. [3] 

was not exactly conserved. 

In the other subsections, the predictions obtained with the new models of the electromagnetic current will be 
compared with data in A=2 and 3 nucleon systems. 



A. Test of the wave functions 



In order to test the PHH wave functions, we have performed a series of calculations of the pd capture reaction 
at Tc.m. = 2 MeV with a Hamiltonian including the Argonne Uig (AV18) two-nucleon [8] and Urbana IX (UIX) 

three- nucleon [9] interactions (the AV18/UIX Hamiltonian model). The model for the electromagnetic current chosen 
for this test is the new "full" one, i.e. including the one-body, the new-ME two-body and the ME three-body terms. 
More precisely. 



i i<j cyclic ijk 

where j^^'^'™{q) and jjfki{q) are given in Eqs. (2.65) and (3.7), respectively. The RMEs are computed from the 
matrix elements 



= i^ftrxi^) -iHcimif"') , (5-2) 

where ex{ci), A = ±1, are the spherical components of the photon polarization vector [2]. 

Some of the most relevant RMEs for pd capture at this energy are the electric dipoles £1 induced by transitions 
between pd. states in relative orbital angular momentum quantum number L = 1,3 and the '"'He state. Here S — 
1/2,3/2 are the channel spin quantum numbers obtained by coupling the spins of the proton and deuteron, and 
J = L + S (the notation and definition used for the RMEs are those of Ref. [3]). The calculated RMEs are listed in 
Table I. In the different calculations, we varied the number of channels included in the bound and scattering PHH 
wave functions, namely the value A^c in the first sum of Eq. (4.5). The channels are ordered for increasing values of 
la + Loc, la and La being the orbital angular momentum of the pair and of the third nucleon with respect to the pair, 
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TABLE I: RMEs ( x 10^) for pd radiative capture at Tc.m. = 2 MeV obtained with the AV18/UIX Hamiltonian and the new-ME 
two- and three-body currents. For the exact definition of the RMEs, see Eq. (4.26) of Rcf . [3]. In the table, A^c(b.s.) (A^c(s.s.)) 
indicates the number of channels included in the expansion of the bound state (scattering state of given J and parity tt). 
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respectively. In this analysis, the values for Nk{o) were taken large enough to have full convergence with respect to 
the order K of Jacobi polynomials. 

First, consider the effect of the truncation of the PHH expansion in the description of the bound-state. The 
calculated binding energy for ^He with the AV18/UIX potential is 7.725 (7.741) MeV after the inclusion of = 12 
(18) channels in the PHH expansion. As can be seen by inspecting the two columns corresponding to the cases 
A'^c(b.s.) = 12 and 18, the changes in the values of the RMEs is at most 2%. We have checked that the inclusion of 
additional channels in the bound state wave function produces tiny changes in the binding energy (less than 10 keV) 
and negligible changes in all the RMEs. 

Next, consider the convergence with respect to A^c(s.s.), the number of channels in the scattering wave functions. 
In general, the dependence of the calculated RMEs on Nc{s.s.) is weak. The only exceptions are the RMEs due to the 
inhibited Ei transitions proceeding through the spin-channel S = 3/2 states, which require the inclusion of a fairly 
large number of channels. In general, the convergence can be checked by looking at Nd elastic scattering phase-shifts 
^LSJ obtained with the given value of Nc{s.s.). Let us consider, for example, the transitions to the J = \ pd 
scattering state. The clastic phase shift S^^"' for the state L = 1, 5 = 1/2 and J = 1/2 were found to have the values 
i5^5 5 — —7.393°, —7.366° and —7.365° for Nc{s.s.) = 10, 14 and 18 channels, respectively. The corresponding changes 
in are given in the first row of Table I and are very small. On the other hand, for the L = 1, S ^ 3/2 and 

J = 1/2 state, the elastic phase-shift turns out to be S'^i^ = 21.666°, 22.318° and 22.319° for A^c(s.s.) = 10, 14 and 18 
channels, respectively. The clastic phase-shift here has a fairly large change passing from Nc = 10 to Nc = 14, due to 

_l31 

the appearance of important channels in the PHH expansion. The corresponding change in S-^^^ is very significant, 
as can be seen by inspecting the second row of Table I. However, by adding more channels, the elastic phase-shift 
and the corresponding capture RMEs show only tiny changes. A similar check has been performed for all the other 
scattering states included in the calculation. 

We have verified by direct calculation that differences between the observables obtained by using the RMEs com- 
puted with the two largest values of Nc are completely negligible. Note that in the present paper the bound and 
scattering wave functions have been obtained on a more extended grid and with more PHH components than in 
previous publications [2, 3]. However, this better accuracy in the wave functions has produced negligible changes in 
the nd and pd capture observables for Tcm < 2 MeV, which were the focus of Refs. [2, 3]. 

For the range of energies considered here, the most important Nd scattering waves are those with J'^ = ^ , | 
and I . For these scattering states a fairly large number of channels has to be included in the PHH expansion of the 
"core" wave function ^'g''^. Scattering states with higher values of J'^ give very small contributions. As mentioned 
before, we have retained the full PHH expansion in the states up to Jmax = 7/2. For larger values of J, the scattering 
wave function has been approximated as in Eq. (4.11). 

B. Test of the two- and three-body current models 

We have calculated the T20 and T21 observables at Tc.m. = 2 MeV using the Argonne vq (AV6) [50] , the Argonne vs 
(AV8) [51] and AV18 two-nucleon interaction. The AV6 interaction is momentum-independent, while the momentum- 
dependence of the AV8 is due only to the spin-orbit operator. The results are shown in Fig. 6. 
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First, let us consider the calculations performed with the AV6 interaction, reported in panels (a) and (b). The 
dotted curves are obtained by including only the one-body current contributions. The dashed curves are obtained 
when the contributions of the old-ME two-body currents of Refs. [2-4], see Eq. (2.64), are added to the one-body 
ones. The solid curves are obtained including instead the two-body current contributions calculated within the MS 
scheme and using the linear path (LP-MS, see Sec. II B). Finally, the dotted-dashed curves are obtained in the long- 
wavelength-approximation (LWA). We observe that in this case there is no significant difference between the old-ME, 
LP-MS, LWA and experimental results. For this potential v^j = 0, and therefore the old-ME and new-ME two-body 
current models coincide. We see that the current derived from the momentum-independent part of the interaction 
using cither the ME or LP-MS schemes are almost equivalent. This is not too surprising, given the small photon 
energy involved in the process, q « 0.035 fm~^, see discussion in Sec. II B. 

Next, let us consider the calculations performed with the AV8 and AV18 interactions, reported in panels (c)-(f). 
Now, the solid curves are obtained using the new-ME scheme, namely the current jf^^''^^^{q) given in Eq. (2.65). 
Note that the old-ME two-body current model results (dashed lines) are in significant disagreement with the LWA 
ones and the experimental data, as can be seen by inspecting panels (c)-(f) of Fig. 6. This is not the case for the 
new-ME results (those obtained with the LP-MS model have not been reported, since they are practically coincident 
with the solid lines). This indicates that the current operator j^^{q]v^) used in Refs. [2-4] and in earlier studies, 
which does not satisfy exactly the CCR with the momentum-dependent terms of the two-nucleon interaction, contain 
spurious contributions. 

This conclusion is supported by another observation. Using the bound and scattering wave hmctions derived from 
the AV18 interaction, but taking into account jj^^{q;v^) only, a good description of the observables T20 and T21 is 
still obtained (see, for example, Fig. 24 of Ref. [24]). The inclusion of the old-model current jif'^ {q;vP) produces 
quite large effects on the RMEs and spoils such an agreement with the data [3] . On the other hand, the contribution 
of the momentum-dependent part of the interaction is noticeably smaller than that one produced by v^j , as can be 
seen, for example, in studies of the binding energies of the light nuclei [52]. Therefore, one can reasonably expect 
that also \jij{q; v^)] <C \jij{q', «")!• The current j[j^iq] v^), constructed in order to properly satisfy the CCR with vfj, 
gives correctly a small contribution to the RMEs, and now the T20 and T21 at Tc.m. = 2 MeV are well reproduced. 
The same happens at higher energies, as will be shown below. 

In order to verify that the agreement found in Ref. [3] for other observables is not spoiled, the differential cross 
section, proton vector analyzing power and the four deuteron tensor analyzing powers for pd capture at Tc.m.= 2 and 
3.33 MeV, calculated with the AVIS two-nucleon interaction, are compared with the experimental data of Refs. [20, 49] 
in Figs. 7 and 8. In the figures, the dotted curves are obtained with only one-body current contributions, the dashed 
and dotted-dashed curves are obtained using the one-body plus "model- independent" (MI) contributions jf^^'°^'^{q) 
and Jj*^^'"'^™(g), respectively, and the solid curve is obtained when, in addition to jf^^'^™ {q)^ the "model-dependent" 
(MD) contributions {jij^{q)), due to the pnj and luttj transition currents and to the current associated with the 
excitation of one intermediate A resonance, are retained ("full" model). The contributions from the isospin-symmetry- 
breaking operators described in Sec. II C are also included, but have been found to be completely negligible. 

With the new model for the nuclear current operator, there is an overall good agreement between experimental 
results and theoretical predictions, except for the iTn observable at small cm. angles. Furthermore, comparing the 
solid and dotted-dashed lines, we conclude that the MD contributions are typically very small, the only exception 
being those for iTn. As will be shown below, this observable is also influenced by three-body current contributions. 
The improved description of the measured T20 and T21 observables, discussed earlier, is evident. 

We now turn our attention to the three-body current. In Fig. 9 the tensor spin observables T20 and T21 for pd 
radiative capture at Tc.m. = 2 and 3.33 MeV are calculated using the wave functions from the AV18/UIX Hamiltonian 
model. The dotted curves correspond to calculations with one- and new-ME two-body currents only. The dashed 
curves have been obtained by including, in addition, the three-body current of Ref. [4], obtained within the old-TCO 
approach. Finally, the solid curves correspond to calculations with one-, new-ME two-body and new-ME three-body 
current jfj^{q) of Eqs. (3.22) and (3.7), obtained within the ME scheme. The results obtained with the three-body 

current operator calculated within the LP-MS scheme {j^j^{q) of Eq. (3.18)) are not shown, because they coincide 
with those obtained with the ME method. Finally, the dotted-dashed curves are the LWA results. Inspection of 
Fig. 9 indicates that, if we use the wave functions obtained from a Hamiltonian including a TNI but disregard the 
corresponding three-body current (dotted curves), there is a significant disagreement with the data. The use of 
the old-TCO three-body current improves partially the description of the data, but only including the new-ME (or, 
equivalently, the LP-MS) three-body current leads to a satisfactory agreement with the experimental data and LWA 
results, especially for the T21 observable. 

Finally, the differential cross section and the spin polarization observables for pd capture at Tc,m.= 2 and 3.33 
MeV, calculated with the AV18, AV18/UIX and AV18/TM Hamiltonian models are compared with the experimental 
data of Refs. [20, 49] in Figs. 10 and 11. The dashed lines are the AV18 results obtained using the corresponding 
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TABLE II: Isoscalar and isovector combinations /xs and of the ^He and magnetic moments, in nuclear magnetons, 
compared with experimental data. The results labeled "lb" are obtained with single-nucloon currents only, those labeled "full- 
new" retain in addition two- and three-body currents in the new model summarized in Sees. II D and IIIC. Also listed are the 
results obtained with the old-ME two-body and old-TCO three-body currents of Ref. [4] ("full-old"). The experimental data 
are from Ref. [53]. 





lb 


0.407 2.165 


Full-new 


0.414 2.539 


Full-old 


0.442 2.557 


Expt. 


0.426 2.553 



wave functions and including, in addition to the one-body current operator, the new- ME current jf^ '"""^{q) and the 

MD current jf-'^{q). The thin solid curves arc obtained with the Haniiltonian including the AV18 and the Tucson- 
Melbourne (TM) [39] TNI (AV18/TM model). In this case, the current includes the one-body, the new-ME and MD 
two-body and the new-ME three-body currents (constructed to satisfy the CCR with the AV18/TM Hamiltonian). 
The thick solid lines arc finally obtained with the AV18/UIX Hamiltonian, and include the corresponding set of one-, 
two-, and three-body currents in the new-ME scheme. As the figures suggest, there are no significant differences 
between the AV18, AV18/TM and AV18/UIX results, except for some tiny effects in the iTn observable. 

C. Magnetic structure of A—S nuclei 

The isoscalar and isovector combinations of the magnetic moments and form factors of ^He and are given in 
Tabic II and Figs. 12 and 13, respectively. The nuclear wave functions have been calculated using the the AV18/UIX 
Hamiltonian model. The results labeled "lb" are obtained retaining only the one-body current operator, those 
labeled "full-new" are obtained including, in addition, the new-ME two-body current contributions and the three- 
body current contributions calculated in the ME scheme. Also listed, arc the results obtained with the old-ME 
two-body and old-TCO three-body currents, as in Ref. [4], labeled "full-old". These last results are slightly different 
from those reported in Ref. [4] , due to the present use of more accurate trinucleon wave functions. The experimental 
data are from Refs. [53 62]. 

Note that: (i) the "full-old" and "full- new" results for the isoscalar and isovector magnetic moments differ by less 
than 1 % and arc very close to the experimental data; (ii) the experimental results for the isovector magnetic form 
factor are fairly well reproduced for momentum transfer < 3.5 fm~^, and the "full-new" curve is slightly closer to 
the experimental data in the region > 4 fm~^ than the "full-old" curve; (iii) the "full-new" curve for the isoscalar 
magnetic form factor is closer to the experimental data than the "full-old" curve in the region < 4 fm~^; (iv) the 
"full- new" curves for the isoscalar and isovector form factors are in disagreement with the data for > 4-4.5 fm~^, 
and the discrepancy between theory and experiment in this region remains unresolved. However, the experimental 
data for the isoscalar magnetic form factor have large errors at high values. 

D. A— 2 radiative capture reaction and deuteron photodisintegration 

The calculated values for the ^H(n,7)^H cross section at thermal neutron energies with the old- and new-ME models 
of the current are listed in Table III. The AV18 two-nucleon interaction is used. 

The small difference between the results obtained with the old- and new-ME models is due to differences in the 
isovector structure of the two-body currents from the momentum dependent terms of the AVIS. The result with the 
new-ME model happens to be in perfect agreement with the experimental value reported in Ref. [63]. 

The deuteron photodisintegration cross sections up to 20 MeV photon energies obtained in impulse approximation 
and with the full-old and full- new ME current models are shown in Fig. 14, along with the experimental data [64-70] . 
Also shown in Fig. 14 are the predictions in which the dominant Ex transitions connecting the deuteron and np 
triplet P-waves are calculated using the Siegert form for the E\ operator, valid in the LWA limit. Again, differences 
between the results obtained with the new- and old-ME current models is to be attributed mostly to differences in 
the isovector currents originating from the momentum dependence of the AV18. Indeed, these terms ensure that 
the new-ME current is exactly conserved, and make the corresponding results essentially identical to the Siegert 
predictions. 
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TABLE III: Total cross section in mb for the np radiative capture, calculated using the AV18 two-nucleon interaction. The 
results labeled "lb" are obtained with single-nucleon currents only, those labeled "lb+2b-MI(old-ME)" and "lb+2b-MI(new- 
ME)" retain in addition model-independent two-body currents in the old-ME and new-ME model summarized in Sees. II D. The 

results labeled "full-old" and "full-new" are obtained by adding the contributions of the model-dependent two-body currents 
to the "lb-l-2b-MI(old-ME)" and "lb-|-2b-MI(new-ME)" results, respectively. The experimental value is from Ref. [63]. 

a{mh) 

lb 304.6 

lb-H2b-MI(oId-ME) 326.1 
lb-F2b-MI(new-ME) 324.7 
Pull-old 334.2 

Full-new 332.7 

Expt. 332.6 ± 0.7 



TABLE IV: Cross section (in mb) and photon polarization parameter Rc of the reaction ^U{n, 7)'^H for the AV18/UIX potential 
model at thermal neutron energy. The results labeled "lb" arc obtained with siugle-uucleon currents only, those labeled 
"lb-t-2b-MI(old-ME)" and "lb-|-2b-MI(new-ME)" retain in addition MI two-body currents in the old-ME and new-ME scheme, 
respectively. The results labeled ". . .-|-2b-MD" are obtained by adding the model-dependent two-body currents to the single- 
nucleon and the MI new-ME two-body currents. Finally, the results labeled "full-new" are obtained by including to the latter 
model also the the contribution of the ME three-body currents, defined in Sec. IIIC. The experimental values for ar and Rc 
are from Refs. [71, 72], respectively. 



Current component 


(7t 


Rc 


lb 


0.227 


-0.061 


lb-|-2b-MI(old-ME) 


0.462 


-0.446 


lb-|-2b-MI(ncw-ME) 


0.418 


-0.429 


. . .-h2b-MD 


0.523 


-0.469 


Full-new 


0.556 


-0.476 


Expt. 


0.508 ± 0.015 -0.42 ± 0.03 



E. A—S radiative capture reactions 

We report here the results for the radiative capture reactions ^H(n,7)'^H and ^H(p, 7)'^He, obtained with the 
AV18/UIX Hamiltonian model. 

1. The '^H(n,^)^H radiative capture reaction 

At thermal energies the nd capture reaction proceeds through S'-wave capture predominantly via magnetic dipole 
transitions from the initial doublet J=l/2 and quartet J=3/2 nd scattering states. In addition, there is a small 
contribution due to an electric quadrupole transition from the initial quartet state. 

The results for the thermal energy cross section and photon polarization parameter are presented in Table IV, 
along with the experimental data [71, 72]. As can be seen by inspection of the table, the cross section calculated with 
single-nucleon currents is approximately a factor of 2 smaller than the measured value. A previous calculation [2] gave 
a-T(lb) = 0.223 mb, very close to the results presented in the first row of Table IV. Inclusion of the model-independent 
new-ME two-body currents leads to a value of (Tt 10% smaller than obtained earlier with the MI old-ME currents of 
Ref. [2]. By adding the model-dependent two-body current, an estimate of ctt = 0.523 mb is obtained. This value 
is to be compared with the corresponding result ctt = 0.558 mb obtained in Ref. [2]. The use of the present MI 
two-body current operators therefore leads to an estimate closer to the experimental datum ctt = 0.508 ±015 mb [71]. 
However, the addition of the three-body currents, which give a rather sizable contribution as can be seen from the 
row labeled "full-new" in Table IV, brings the total cross section to ctt = 0.556 mb. The 9% slight overprediction is 
presumably due to the model-dependent currents associated with the A excitations. Fortunately, at Tc.m. > 1 MeV, 
this MD current gives a negligible contribution to the cross section and the other polarization observables, as already 
shown in Figs. 7 and 8. 

The photon polarization parameter is very sensitive to two-body currents (for its definition in terms of RMEs, 
see Ref. [2]). For example, for the AV18/UIX Hamiltonian, their inclusion produces roughly a six- fold increase, in 
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absolute value, of the single-nuclcon prediction. Also in this case, we find a 13% overprediction (in absolute value) 
of this parameter. The small reduction of \Rc\, found when the new-ME model for the two-body current is used, is 
compensated by the inclusion of the three-body currents. 

At higher energies, there exist several measurements of the unpolarized differential cross section for both the 
radiative capture process ^H(n, 7)^H [73] and for the "time-reversed" process ^H(7,n)^H [82-86]. In the cm. system 
and at low energies, the unpolarized cross sections are related by the principle of detailed balance 

where q and p are the 7 and the relative nd momenta, respectively. In Figure 15, we compare our predictions for 
{d(Jca.pt/di^)c.m. with the experimental results of Ref. [73], which is the only direct measurement of the ^H(n,7)'^H 
differential cross section. The dashed curves represent the results obtained with the inclusion of the single-nucleon 
current only, the dotted-dashed curves are obtained with one- and new-ME two-body contributions (both MI and 
MD), and the solid curves represent the "full- new" result, obtained including, in addition to the one-body and 
above-mentioned two-body currents, also the ME three-body current contributions. At these energies the process is 
dominated by Ei transitions between the nd P-wave states and the '^H ground state, as can be inferred from the 
bell-shape of the curves. The slight distortion of the peak is due to non-negligible contributions from RMEs, 
coming, in particular, from the J = 3/2,5/2 states with S = 3/2. Our "full-new" calculation reproduces quite well 
the experimental data, with some differences at large angles. As will be shown below, this has some consequences for 
the so-called fore-aft asymmetry, discussed in Sec. VE3. 

In Fig. 16, we compare our "full-new" results with those obtained in Refs. [25, 26]. In Ref. [25], the same potential 
model (AV18/UIX) as in present work has been used, but the authors adopt a slightly different current model (they 
do not consider jfj^{q;vP), jfj^{q) and the three-body current). In Ref. [26], the exchange currents are taken into 
account using Siegert's theorem, and a different two-body potential model has been used (Bonn A [87]), without 
any inclusion of three-nucleon forces. As can be seen by inspecting Fig. 16, the three theoretical calculations are 
practically the same, with some differences with results of Ref. [26] at Tc.m. = 7.20 MeV. This difference is likely due 
to the use of a different potential model, which slightly underestimates the '^H binding energy. 

In addition to unpolarized cross sections, there are also a few analyzing-power angular distribution data [73], but 
they have large error bars, and therefore we have decided not to perform a comparison for this observable. 



2. The ^ H(p,'-^f' He radiative capture reaction 

In Fig. 17 we present the results for the astrophysical S factor of the ^H(p, 7)^He radiative capture reaction at 
thermal energies. This quantity is defined as 



S{T,,m) = T,,^. cJT(Tc.m.) e^^"/""' , (5.4) 

where Tc.m. is the pd cm. kinetic energy, axiTcm.) is the total cross section, a is the fine structure constant, and Vrei 
is the pd relative velocity. The experimental data are from Refs. [16, 88-90]. The solid curve represents the "full-new" 
result, obtained including, in addition to the one-body currents, also the new-ME two-body current contributions and 
the three-body current contributions calculated in the ME scheme. The dashed curve represents the result obtained 
with the inclusion of the single-nucleon current only. Here, no significant difference has been seen between the results 
obtained with the present model for the nuclear current operator and the "old" one of Refs. [2, 3]. The agreement 
between the theoretical predictions and the experimental data, especially the very recent LUNA data [16], is excellent. 
In particular, the calculated S factor at zero energy is 0.219 eV b [3], in very nice agreement with the LUNA result 
of 0.216±0.010eV b. 

There exist several measurements of ^H(p,7),''Hc obscrvables between Tc.m. = 2 and 20 MeV. In Figs. 18-23, we 
compare the predictions obtained with our new model of the current with a selected set of observables. In all these 
figures, the dashed lines are obtained with only one-body contributions, the dotted-dashed ones are obtained with 
one- and new-ME two-body contributions (MI+MD), the solid curves are the "full" results with, in addition, also 
three-body contributions, obtained in the ME scheme. For completeness, the predicted angular distributions of the 
differential cross section a/ao, proton vector analyzing power Ay, and deuteron vector and tensor analyzing powers 
iTii, T20, T21, T22, at Tc.m. = 2 and 3.33 MeV are again given in Figs. 18 and 19. Note that these two cm. energies are 
just below and above the deuteron breakup threshold. We can conclude that: (i) an overall nice description for all the 
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obscrvablcs has been obtained, with the only exception of the iTn deutcron polarization observable at small angles; 
(ii) some small three-body cmrents effects are noticeable, especially in the T20 and T21 deuteron tensor observables. 

The predicted angular distributions of the deuteron vector and tensor analyzing powers Ay{d), Axx, ^yy, A^z, at 
Tc.m. = 5-83 are given in Fig. 20. The experimental data are from Ref. [91]. Comments similar to those above can be 
made in this case too. 

The differential cross section da/dfl and the deuteron vector- and tensor- analyzing powers Ay{d) and Ayy arc given 
in Figs. 21, 22 and 23 for four different cm. energies. The differential cross section is nicely reproduced by theory at 
Tc.m. =6.60 and 9.86 MeV, while some discrepancies are present at Tc.m. =16.00 and 18.66 MeV. However, it should 
be pointed out that these data sets are quite old, and new experimental studies of this process in this energy range 
would be very useful. The Ay{d) observables are poorly reproduced at small angles for Tc.m. =5.83 and 9.66 MeV. 
However, at Tc.m. =15.00 MeV, the discrepancy between theory and experiment seems to disappear. It would be 
interesting to continue this comparison at higher values of Tc.m.- The Ayy observables arc nicely reproduced in the 
whole range of Tc.m.- Some small discrepancies are present at small angles for Tc.m. =15.00 MeV. It is important to 
note, however, that the Ay{d) and Ayy observables are obtained dividing for the differential cross section, which is 
close to zero at small and large values of the cm. angle. In view of this, the agreement between theory and experiment 
for the Ay{d) and Ayy observables should be considered satisfactory in the whole range of cm. angles. 

Finally, in Fig. 24 we study the importance of including the Coulomb interaction in the bound- and scattering- 
state wave functions. In fact, the Tc.m. =3.33 MeV differential cross section and vector- and tensor-analyzing powers 
are calculated using the AV18 nuclear Hamiltonian and one-body plus new-ME two-body currents. The Coulomb 
interaction is included both in the bound- and scattering-state wave functions (thick-solid lines), in the boimd- but 
not in the scattering-state wave functions (thin-solid lines), and neither in the bound- nor in the scattering-state wave 
functions (dashed lines). The Coulomb interaction plays a small but significant role, particularly in the differential 
cross section. 



3. The fore-aft asymmetry 



A quantity of particular interest, both for historical reasons and for purposes of comparison with data, is the 
so-called fore-aft asymmetry in the angular distribution of the cross section. This quantity is defined according to 



a(54.7°)-a(125.3°) 
^ ~ C7(54.7°) + <7(125.3°) ' ^ ' 

where (T{6c.m.) ~ {da ca.pt /d^) cm. and Oc.m. is the ^-N scattering angle. The ^H(p, 7)^He and ^H(n,7)'^H asymmetries 
resulting from our calculation and from that of Ref. [26] arc compared with existing data in Fig. 25. 

The pd asymmetry is well reproduced by the calculation (it is a consequence of the good agreement between the 
theoretical and experimental cross section angular distributions, shown in Fig. 21). On the other hand, there are large 
differences between the theoretical and experimental nd asymmetries. These differences could be due to problems 
in the analysis of the data to extract the experimental values of a^. For example, the experimental asymmetry 
comes out mainly from the measurement taken at the largest angle (see Fig. 15). Note that for other values of Oc.m.i 
good agreement is obtained between theory and experiment. There are also inconsistencies between the asymmetries 
given in Rcfs. [73, 84], therefore it is likely that these discrepancies are due to experimental problems. However, if 
experimentally confirmed, this problem could be of relevance, since the present calculations seem unable to predict 
\as{nd)\ > 0.1. 



VI. SUMMARY AND CONCLUSIONS 



Wc have investigated two different approaches for constructing conserved two- and three-body electromagnetic 
currents: one is based on meson-exchange mechanisms, while the other uses minimal substitution in the explicit and 
implicit-via the isospin exchange operator-momentum dependence of the two- and three- nucleon interactions. In the 
meson exchange model developed by Riska and collaborators [10, 12, 32, 33] and used in earlier studies [2, 3, 11 13], 
some of the terms associated with the isospin- and momentum-dependent components of the two-nucleon interaction 
were ignored, since, due to their short range character, they were expected to give negligible contributions. The 
resulting currents, however, were not strictly conserved. This limitation is removed in the present work. We have 
also shown how the two-body currents obtained in the meson-exchange framework (more precisely, their longitudinal 
part) can be derived in the minimal-substitution scheme, originally developed by Sachs [34], by selecting a specific 
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path in the spacc-cxchangc operator. Lastly, wc have constructed a reahstic model for the three-body electromagnetic 
current satisfying the current conservation relation with the Urbana or Tucson- Melbourne three-nucleon interactions. 

A variety of observables have been calculated to test the present model of nuclear current operator. In particular, 
for the A=3 nuclear systems, cross sections as well as polarization observables have been calculated and compared 
with the corresponding experimental results in the energy range 0-20 MeV. 

In general, the contributions due to the new two- and three-body currents — from the momentum dependence of 
the two-nucleon interaction and from the three-nucleon interaction, respectively — are found to be numerically small. 
However, they resolve the discrepancies between theory and experiment obtained in earlier studies [3] for some of the 
polarization parameters measured in pd radiative capture, specifically the tensor polarizations T20 and T21. These 
contributions also reduce the overprediction of the nd radiative capture cross section at thermal neutron energy from 
the 15 % obtained in Ref. [2] to the current 9 %. 

In conclusion, (i) the predictions for the np radiative capture and low-energy deuteron photodisintegration [38] , and 
for the magnetic form factors of '^He and "'H, have remained essentially unchanged from those reported in previous 
studies [4]; (ii) a satisfactory agreement between theory and experiment for pd radiative capture observables above 
deuteron breakup threshold up to 20 MeV has been found, in particular for the tensor observables. Some discrepancies, 
however, persist in the vector polarization observables at forward angles. 
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APPENDIX A 



In this appendix we show how the meson-exchange (ME) two-body currents (rather, their longitudinal components) 
from the static part of the two-nucleon interaction can be derived in the minimal substitution (MS) scheme by selecting 
a specific path in the space-exchange operator. To this end, we first write 



V2,ij Ti ■ Tj = ^(«ps,,j + + ^^ ' , (Al) 

a 

Vps,ij = -fpsJ'^r ■ ■ V,) Ya{r) , (A2) 



V.,, = -fv.M X v.) • (cT, X V,) y„(r) , (A3) 

VVS,^, = fvS,aYa{r) , (A4) 

W = , (A5) 

where fps,a, fv,a, fvs,a and rua are appropriate parameters, see Eqs. (2.26)-(2.28). The two-body current which 
satisfies the current conservation relation with V2,ij is then written in MS scheme as 



3tj {q; PS) + jfj iq;V)+ j". (g; VS) 



(A6) 



where each term of the sum is given in Eq. (2.47), with V2,ij replaced by the corresponding PS, V or VS potential. 



First, we consider the simple spin- independent current jfj(^q',VSy. by choosing 'yij 
jfj {q; VS) can be written as 



-7'-^ and using Eq. (2.50), 



jtj{q;VS) = f^s,aya{r)3Uq) ' 

jtM) = Gl{ql){nXT^)^ [ dse"^-^ 



(A7) 



where the path is selected to be 
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Q 



La{x) - ma 



1 1 

2 - - 2 ' 



Rij = ^i^i + rj) , s(x = -^j=ri, s(x = +^]=rj, 



La{x) 



(A8) 



It is straightforward to verify that Ylajiji^'- ^'^) reduces to jij{q; VS) of Eq. (2.34), obtained in the ME scheme. In 
this case, both the longitudinal and transverse components of jij{q; VS) are exactly reproduced. 
Next, using Eq. (2.47), the PS current reads 



(A9) 



which can be rewritten as 



i«.(g)(ai-V,)(<T,-V,)F„(r) 



= (<T,-V,)(a,-V,-)j5(Q)F„(r) 

- [{<Ti-Vi)3t,{qmcT,-V,)Ya{r)] 

- [{rryV,)jt,{q)][{cT,-V,)Ya{r)] 

- i;(r)(^.-V.)(^,-V,)j«(q) . 



(AlO) 



The a-meson in-flight current is again exactly reproduced (first line of Eq. (AlO)). However, only the longitudinal 
components of the contact terms are reproduced by Eq. (AlO) (second and third lines), since, via Eq. (2.49), 



(tTi- Vi)9-j?.(q) = -cTi-qe'"-'-* , 



(All) 
(A12) 



and the last term of Eq. (AlO) vanishes when dotted with q. Therefore, using Eqs. (A10)-(A12), Eq. (A9) can be 
rewritten as 



PS,a 



{ai-Vi){cTj-Vj)jUq)Ya{r) 



+ GUql) in X X,), [e'-J-^rr, (<t, • V,) - e'^-^^a, {a, ■ V,)] r„(r) 

+ additional transverse terms . 

It would be interesting to evaluate the contributions of these additional transverse terms. 
Similar considerations are valid for jfj(q; V), which can be written as 



(A13) 



JtMV) = -/^,,(^[(^.x V,)-(<T, X V,)] j-(g)i; 

^(^. X V.) [{a,xV,)-jUq)]Ya{ 



(r) 



-(^, X V,) [(^.x V0-j.'^(g)]r4r) 

- G^(g^)(TiXT,),[e'<' '-'^.x(^, x V,) 

-e^^-^^cTj X (<7i x Vi)]F„(r) 
+ additional transverse terms . 



(A14) 
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APPENDIX B 



We list here the two-body currents obtained with the MS method from the quadratic-momentum-dependent terms 
of the interaction. To this end, the {L ■ S)^ term is written as 



{cTi ■ L) {o-j ■ L) + {a-j ■ L) {cTi ■ L) 



and the potential functions associated with the L ■ S and L"^ operators are re-defined accordingly, for example 



Vh{r) = Vbir) - -Vbbir) , 



Vg{r) = Vg{r) + -Vbb{r) 



(Bl) 

(B2) 
(B3) 



and similarly for Vbrir) and Vqrir). The spin-orbit currents are then those given in Eqs. (2.57) and (2.60) with 
Vb{r)-,VbT{r) replaced by Vb{r),VbT{r). 
For the L'^ terms, given by 



Vq = li^ [Vq{r) + Vqa{r) (Ti ■ (Tj + Vqr{r) Ti ■ Tj + Vqarir) CTi ■ (TjTi ■ Tj] , 

using the linear path of Eq. (2.51), we find 



(B4) 



jij{q;LL) = 

jf/'iq^LLT) = 



Vq{r) +Vqa{r) CTi ' CTj 
Vgrir) + Vqarir) CTi ■ (Tj 



ir — r X L ) P_ -h - (r X q) X r P+ 



ir — r X L ) i?_ + -(r x q) x ri?+ 



+ iGl{ql){nXTi),{L\rUi{q)} 
where r = Vij, {■ ■ ■} denotes the anticommutator, and fij{q) is defined in Eq. (2.53). We have also defined 



(B5) 



(B6) 



with Ci, Cj, r]i and r]j listed in Sec. II. 

For the quadratic spin-orbit terms, given by 



m = -^Vbb{r){cri ■ L,<Tj-L} 

+ \vbbT{r){{cri ■ L, tTj ■ L} ,Ti- tj} , 



again using the linear path, we find 



(B7) 
(B8) 



(B9) 



jijiq;bb) = -Vbbir) \ P- , [(cr, x r) [ct^ ■ L) + [a-j x r) {(Ti ■ L)] 



jt/{q;bbT) = -Vbbrir) 



{(Tj xr) [q- {cTi X r)] -|- (ct^ xr)[q- {(Tj x r)] 



R. 



{aj X r) {ai ■ L) + {(Ti x r) (cr^ • L) + i{(Ti ■ (Tj)r - - [ai {a^ ■ r) + crj {ai ■ r)] 



+ ^G^{<ll){n X ■L,(Tj L}, rfij{q)^ 



(BIO) 



(Bll) 



APPENDIX C 

Using Eq. (2.33), the configuration-space expressions for the exchange currents jlj{ki,kj; B) of Sec. Ill A, B 
or V, are given by: 



e'"-^* gps{r)ai (S,- • f) + e'""-^' gps{r)Sj {(n ■ r) 



+ e 



iq-R 



.Gps,2ir 



a-i {Sj ■ r) + Sj (tTj • f) + f (o-j • Sj) 
(Tj (Sj • q) - 1 '-^ Sj {cT, ■ q) 



+1 



_.Gps^ . . .^(g^. . ^ . Gps^ . , ^^^g^, , „^ 



-Gps,6(r) r (o-i • g)(S,- • 9) - f {en ■ r)(S^ • r) 



e^""-* gv{r)ai x (S,- x r) + e^""-^ 9vir)Sj x (o-^ x r) 



( (Sj xr)xai + {ai x f) x Sj + 2f (tri • Sj) 



+1 



r 



r 

■ GvAir) 



(Sj X q) X ai - 1 — ((Tj X g) X Sj 

f {ai X f) ■ (Sj X q) + i— ^^^^ f (o-j x q) ■ (Sj x r) 



-Gy,6(r) r {ai x q) ■ (Sj x q) ^^^^^^ r (o-j x r) ■ (Sj x r) 



- r 



iq-R 



Gv.2{r) (Sj X q) (t, • (qr X f) + Gy,3(r) (cr, x q) Sj • (q x r) 



+1 



.Gy,4(r) + Gy,5(r) 



(Sj xr) ai-{fxq)- {ai x r) Sj ■ {r x q) 



+1 — j Sj X (<Ti X q) - (Ti X (Sj X q) 



where r = — Tj, r = r/r, il = \{ri + rj), and 



6'ps(0 = - — 



1 

3r 
1 
3r 

1/2 



dr'r"vii{r') - [ 

J r J r 



dr' 



EB{x;r) - r — EB{x;r) 
ar 



GB,i{r) = I daje-'^'" 
-1/2 

"1/2 1 

GB,2(r) = / dxe-'-''-'^{-+x)EB{x;r) 

-1/2 ^ 

»l/2 

Gs,3 

■ -1/2 



d 



5(r) = ['^' dxe-'-'' '^{l--x)EB{x;r) 
J-1/2 ^ 
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Gs,4(r) 



f^'^ ■ 1 
j dxe-'^'i-^{-+x) 



EB{x;r) 



/2 

1/2 . 

dxe-'^«-''(- - a;) 

-1/2 2 
1/2 -, 

■1/2 4 



r — EB{x;r) 
ar 



EB{x;r) - r — EB{x;r) 



1/2 



1/2 



d d^ 
3EB{x;r) - 3r — EB{x;r) + EB{x;r) 



The functions EB{x;r) are defined as 



(C8) 
(C9) 
(CIO) 
(Cll) 



EB{x;r) = f^^e-^^(-) 



a=l 



K{x) = ^lml + ^{l-Ax^). 



(C12) 



(C13) 



The coefficients gs a are obtained by fitting the functions Vg {k) of Sec. Ill A, B = PS, V, with J2a=i 9b a/(^^ 
ml) [12]. 
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FIG. 1: Integration paths from position i to position j and vice versa for the integral of Eq. (2.40). 
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FIG. 2: Feynman diagram representation of the three-body currents associated with PS- and V^-meson exchanges. Thin, thick, 
dashed and wavy lines denote nucleons, A-isobars, a-mesons, and photons, respectively. 





FIG. 3: Integration paths chosen to derive the three-body current operator of Eqs. (3.9) and (3.14). 




FIG. 4: Final position of particles i, j, and k after that the product of the space exchange operators PjkPij is applied. 
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e [deg] [deg] 

cm. cm."- 

FIG. 5: Differential cross section, proton vector-analyzing power Ay, and four deuteron tensor polarization observables for pd 
elastic scattering at Tc.m. =6.66 MeV as function of the cm. scattering angle. The dashed and solid curves are obtained with 
the AV18 and AV18/UIX interaction models, respectively. The experimental differential cross section and Ay data are from 
Ref. [46], while the experimental tensor-analyzing power data are from Ref. [47]. 
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FIG. 6: Deuteron tensor polarization observables T20 and T21 for -pd radiative capture at T^.m. =2 MeV as function of the 
cm. 7-p scattering angle, obtained with the AV6, AV8 and AV18 Hamiltonian models. Dotted, dashed and solid lines are 

obtained with only one-body current, with one- and old-ME two-body currents and with one- and LP-MS two-body currents, 
respectively. The results obtained in long- wavelength approximation are also shown (dotted-dashed lines). The experimental 
data are from Ref. [49]. 
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FIG. 7: Differential cross section, proton vector analyzing power, and the four deuteron tensor polarization observables for pd 
radiative capture at Tc.m. =2 MeV as function of the cm. 7-p scattering angle, obtained with the AV18 Hamiltoniam model. 
The dotted curves axe obtained with only one-body currents, the dashed and dotted-dashed curves retain, in addition, the 
contributions from the MI two-body operators and j^^'^^^^iq) (see Eqs. (2.64) and (2.65)), respectively. The solid 

curves arc obtained including, in addition to j^j^''^'^^ (q) , also the "model-dependent" two-body contributions. The experimental 
data are from Ref. [49]. In most of the panels, the dotted-dashed and solid curves are indistinguishable. 
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FIG. 8: Same as Fig. 7 at Tc.m. =3.33 MeV. The experimental data are from Ref. [20]. 
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FIG. 9: Deuteron tensor polarization observables T20 and T21 for pd radiative capture at Tc.m. =2 and 3.33 MeV as function 
of the cm. 7-p scattering angle, obtained with the AV18/UIX Hamiltonian model. The dotted, dashed and solid curves 
correspond to the calculation with one- and two-body currents only, with the three-body current obtained within the TOO 

approach, with the thrcc-body current jijif(q) obtained within the ME scheme (see Eqs. (3.22) and (3.7)). The results obtained 
in LWA are also shown (dotted-dashed curves). The experimental data are from Refs. [20, 49]. 
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FIG. 10: Differential cross section, proton vector analyzing power, and the four deuteron tensor polarization observables for 
pd radiative capture at Tc.m. =2 MeV as function of the cm. j-p scattering angle, obtained with the AV18 (dashed lines), 

AV18/TM (thin solid lines) and AV18/UIX (thick solid linos) Hamiltonian models. The model for the nuclear current operator 
include one-, two-, and three-body contributions and they satisfy the CCR with each given Hamiltonian. The experimental 
data are from Ref. [49]. In most of the panels, the three curves are indistinguishable. 
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FIG. 11: Same as Fig. 10 at Tc.m. =3.33 MoV. The dashed and thick sohd hnes are obtained with the AV18 and AV18/UIX 
Hamiltonian models, respectively. The experimental data are from Ref. [20]. 
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FIG. 12: The isoscalar combination of the ''He and "^H magnetic form factors, obtained with single nucleon currents (lb), and 
with the inclusion of two- and three-body currents in the new model summarized in Sees. II D and IIIC (full-new). Also listed 
are the results obtained with the old-ME two-body and old-TCO three-body currents of Ref. [4] (full-old). The experimental 
data are from Ref. [54-62]. 




FIG. 13: Same as Fig. 12, but for the isovector combination of the ^He and "^H magnetic form factors. The experimental data 
are from Ref. [54-62]. 
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FIG. 14: The deuteron photodisintegration cross sections obtained with one-body current (lb) and with, in addition, the old 
(full-old) or new (full- new) models for two-body currents are compared to experimental values. Also shown are the results 
obtained by using the Siegert form for the Ei transition. The experimental data are from Refs. [64-70]. 
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FIG. 15: Differential cross section in /ib/sr for nd radiative capture at Tc.m. =6 MeV and 7.2 MeV as function of the cm. 7-n 
scattering angle, obtained with the AV18/UIX Hamiltoniam model. The dashes curves have been obtained using the single- 
nucleon current only. The dotted-dashed curves have been obtained with the inclusion of the new-ME two-body currents. 
Finally, the solid curves have been obtained by adding the ME three-body current. The experimental data are from Ref. [73]. 
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FIG. 16: Differential cross section in /ib/sr for nd radiative capture at Tc.m. =6 MeV and 7.2 MeV as function of the cm. y-n 
scattering angle, obtained in present work (thick solid curves) and in Ref. [25] (thin solid curves) using the AV18/UIX potential 
model. The two curves are practically indistinguishable. The dashed curves are the results of Ref. [26] obtained with the Bonn 
A potential model [87]. The experimental data are from Ref. [73]. 




FIG. 17: The S factor of the ^H(p, 7)^He reaction, obtained with the AV18/UIX Hamiltonian model and with the one-body 
current only (dashed line), and with the one-, two- and three-body current (solid line) is compared with the experimental 
results from Refs. [16, 88-90]. 
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FIG. 18: Differential cross section, proton vector analyzing power, and the four deuteron tensor polarization observables for 
pd radiative capture at Tc.m. =2 MeV as function of the cm. 7-p scattering angle, obtained with the AV18/UIX Hamiltonian 
model. The dashed, dottcd-dashed and solid lines are obtained with one-body contributions, one- and new-ME two-body 
contributions, and one-, two- and three-body contributions. The experimental data are from Ref. [49]. 
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FIG. 19: Same as Fig. 18 at Tc.m. =3.33 MeV. The experimental data are from Ref. [20]. 
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FIG. 20: Four dcutoron tensor polarization obscrvablcs for pd radiative capture at Ten. =5.83 MeV as function of the cm. 7-p 
scattering angle, obtained with the AV18/UIX Hamiltonian model. The same notation as in Fig. 18 is used for the different 
lines. The experimental data are from Ref. [91]. 
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FIG. 21: Differential cross section in /ib/sr for pd radiative capture up to Tc.m. =18.66 MeV as function of the cm. y-p 
scattering angle, obtained with the AV18/UIX Hamiltonian model. The same notation as in Fig. 18 is used for the different 
lines. The experimental data are from Ref. [75] for Tc.m. =6.60 and 9.86 MeV, and from Ref. [79] for the other cases. 
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FIG. 22: Deuteron vector analyzing power Ay{d) for pd radiative capture up to Tc.m. =18.66 MeV as function of the cm. 7-p 
scattering angle, obtained with the AV18/UIX Hamiltonian model. The same notation as in Fig. 18 is used for the different 
hnes. The experimental data are from Ref. [91] for Tc.m. =5.83 MeV and from Refs. [92, 93] for Tc.m. =9.66 and 15.00 MeV. 
No data at Tc.m. = 18.66 MeV are available in literature up to now. 
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FIG. 23: Deuteron tensor analyzing power for pd radiative capture up to Tc.m. =18.66 MeV as function of the cm. 7-p 
scattering angle, obtained with the AV18/UIX Hamiltonian model. The same notation as in Fig. 18 is used for the different 
Unes. The experimental data are from Ref. [91] for Tc.m. =5.83 MeV, from Refs. [92, 93] for Tc.m. =9.66, and from Refs. [92-94] 
for Tc.m. =15.00 MeV. No data at Tc.m. = 18.66 MeV are available in literature up to now. 
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FIG. 24: Differential cross section, proton vector analyzing power, and the four deuteron tensor polarization observables for pd 
radiative capture at Tc.m. =3.33 MeV as function of the cm. f-p scattering angle, obtained with the AV18 Hamiltonian model 
and the one- plus new-ME two-body nuclear current operators. The dashed, thin-solid and thick-solid lines are obtained with 

no Coulomb interaction in the bound- and scattering-state wave functions, with Coulomb interaction only in the bound-state 
wave function, and with inclusion of the Coulomb interaction both in the bound- and in the scattering-state wave functions, 
respectively. The experimental data are from Ref. [20]. 
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FIG. 25: Fore-aft asymmetry for pd and nd radiative capture as function of the 7 energy, obtained with the AV18/UIX 
Hamiltonian model and the "full- new" current (solid lines). The results of the calculation of Ref. [26] have been also reported 
(dashed lines). The pd experimental data are from Refs. [74-81], and the nd experimental data are from Refs. [73, 82, 84]. 



